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We review current observational constraints on the polarization of the Cosmic Microwave 
Background (CMB), with a particular emphasis on detecting the signature of primordial 
gravitational waves. We present an analytic solution to the Polanarev approximation for 
CMB polarization produced by primordial gravitational waves. This simplifies the calcu- 
lation of the curl, or B-mode power spectrum associated with gravitational waves during 
the epoch of cosmological inflation. We compare our analytic method to existing numeri- 
cal methods and also make predictions for the sensitivity of upcoming CMB polarization 
observations to the inflationary gravitational wave background. We show that upcom- 
ing experiments should be able either detect the relic gravitational wave background or 
completely rule out whole classes of inflationary models. 

Keywords: cosmic microwave background; polarization; measurements 

PACS Nos.: include PACS Nos. 

1. Introduction 

The cosmic microwave background (CMB) is one of the most powerful and precise 
cosmological probes. Because the CMB photons we observe today probe the physics 
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of the early universe during the epoch of hnear gravity, the CMB is often referred 
to as a "snapshot" of the primordial universe. The CMB has the promise to address 
the most fundamental cosmological questions: the geometry and age of the universe, 
the matter-energy content of the universe, the ionization history and the spectrum 
of primordial perturbations. 

This review addresses the theoretical foundations of CMB polarization gener- 
ated by cosmological gravitational waves with particular emphasis given to analytic 
and numerical results which encode its behavior most cogently. Gravitational waves, 
in contrast to adiabatic perturbations (which are the dominant source of CMB tem- 
perature and polariza tion anisotropy) imprint a unique divergence-free pattern of 
polarization on the skjCEl. This pattern is called "B-mode'^^or "cur l-mode'^ polar- 
ization. Although there may be a substantial GWB contribution to CMB tempera- 
ture anisotropy, its effect on the C MB temperature is nearly completely degenerate 
with other cosmological parameter J'^l^l This is unsurprising since the CMB temper- 
ature is a scalar quantity. However, the tensorial nature of CMB polarization permits 
separation of scalar fluctuations from tensor, CWB- gene rated, fluctuations. CMB 
polarization maps can be decomposed into two terma^l^l One term is the gradient 
of a scalar potential and is invariant under parity transforms (often called "E-mode" 
in analogy to the electric field/scalar potential). The second component is the curl 
of a vector potential ("B-mode"). Scalar perturbations have no handedness so the 
primary CMB curl-mode component exists only if there is a GWB. Besides the 
obvious importance of a new method of (indirect) detection of gravitational waves, 
detection of the B-mode signal provides the cleanest, and perhaps only window into 
unique predictions of the inflationary cosmological paradigm. 

Numerical methods for simulating CMB temperature and polarization power 
spectra have revolutionized cosmology . Without such codes parameter estima- 
tion from CMB data sets would be difficult, if not impossible. However, the com- 
plexity of these codes has hidden the underlying physics. For this reason there 
is considerable interest in analytical approaches for calculating CMB polarization 
caused by cosmological gravitational waves (see for example Ref. ^ IHl ^1 

El. 

This review elucidates an analytical approach to the problem using the funda- 
mental physical characteristics of gravitational waves to significantly simplify calcu- 
lations. Our method simplifies the comparison between theoretical predictions and 
future observational results. We compare our results to existing numerical methods 
and summarize current observational results. We conclude with by making predic- 
tions for ground, balloon, and space-based observations in the upcoming decade. 

2. Inflation and Primordial Gravitational Waves 

Inflation "'^'^ is a bold cosmological paradigm which has provided key insights into 
many observations of modern cosmology. Inflation provides an explanation for the 
observed spatial-flatness of the universe, motivates the lack of topological defects 
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Gravitational Waves and the Polarization of the CMB 3 

such as magnetic monopoles, and explains the exquisite isotropy of the cosmic mi- 
crowave background (CMB) while simultaneously providing a mechanism to gen- 
erate the observed fluctuations. However, along with its many successes has come 
increased scrutiny. M any of inflation's key predictions are shared by alternative 
models. While skeptic^^^l are increasingly challenged to attack it, inflation's pro- 
ponents can only claim circumstantial evidence in its favor. Only new discoveries 
will provide the data required to break the deadlock. A conclusive detection of the 
primordial gravitational wave background (GWB) predicted by inflation would be 
"the smoking gun" confirming the inflationary model beyond a reasonable doubt. 
No other known cosmological mechanism mimics the imprint of the GWB on the 
polarization of the CMB. 

One of inflation's successful predictions is its solution to the "horizon problem" 
- the observation that regions of the universe share the same thermodynamic tem- 
perature despite never being in thermal, or even causal, contact. Inflation solves the 
horizon problem via a superluminal expansion of the universe at very early times, 
prior to the ordinary Hubble expansion observed today. This rendered the entire 
observable universe within causal contact initially. This expansion also accounts 

for the (seemingly) finely-tuned spatial flatness of the universe observed by CMB 

1 7 1 S 1 Q 9fl 91 99 

temperature anisotropy experiments . 

This review focuses on inflationary generated gravitational wave, or tensor per- 
turbations. However, inflation also predicts a nearly scale-invariant spectrum of 
energy density, or scalar, perturbations. That any perturbations remain after the 
universe expanded by 60 e-folds is astonishing! Yet, surprisingly, the observed fluc- 
tuation level arises naturally as magnifled quantum fluctuations of the scalar fleld 
that drove inflation (the inflaton) (see |23| for a review of inflationary perturbation 
theory) . Following the inflationary epoch, the size of the residual fluctuations were 
imprinted on the surface of last scattering and are observable in the CMB. The 
observed size of the fluctuations, e.g., and the correlationJ232Sl between the 
CMB's temperature and polarization patterns at super-horizon scales motivates, 
but does not prove, the inflationary paradigm. 

Inflation not only produces these perturbations, but also endows them with 
synchronized initial phases. This is crucial as it allows the perturbations to 
grow without cancellation, as would occur if the initial phases were uncorrelated. 
WMAP an d othe r CMB experiments, in combination with large-scale-structure 
observationJ232Hl probe the scalar perturbation spectrum and are consistent with 
inflation. 

Regrettably, neither flatness nor smoothness are unique to inflation. Both have 
long histories in cosmology. The fl atne ss of the universe was anticipated prior to 
inflation on quasi- anthropic ground J^. A scale-invariant primordial matter/energy 
perturbation spectr um (i.e., d eviations from perfect smoothness) was also predicted 
well before inflatiorP^^^^, though no mechanism to produce the perturbations 
was given in these early works. These perturbations, combined with the universe's 
spatial flatness, increase the plausibility of the inflationary paradigm, since the 
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amount of expansion required for flatness should have also smoothed any initial 
perturbations to zero. 

Two more observations are noteworthy for their consistency with the inflation- 
ary paradigm: the low abundance of rehcts^ (e.g., no magnetic monopoles) and 
th e obser ved Gaussianity of perturbationJ^. Since, as many authors have pointed 
Q^J23SE1 inflation has passed so many observational tests any replacement the- 
ory would need to look very much like inflation. However, while there is abundant 
circumstantial evidence, there is one unique prediction of inflation that cannot be 
mimicked: a primordial gravitati onal wave background; thereby cementing its status 
as "the smoking gun of inflation'^^^. 

Inflation posits a new scalar field (the inflaton) and specifies its action poten- 
tial, and thereby, its dynamics (calculated using the standard methods of scalar 
field theorjE^. All theories of inflation produce a GWB, though some at unob- 
servable levels. All are scalar field theories, incorporating fluctuations via quantum 
perturbation theory. While the identity of the inflaton is u nkno wn, specifying the in- 
flaton potential has immediate observational consequence^^. The most important 
inflationary measure of the GWB is the tensor-to-scalar ratio, r, since it param- 
eterizes the unique prediction of inflation. A detection of r would simultaneously 
reveal both the epoch of inflation and its energy scalc^^. If, as theorists speculate, 
inflation is related to the GUT-scale (grand unified theory) then detection of the 
B-mode signature would probe physics at at the 10^^ GeV scalJ23 



3. CMB Polarization 

The CMB has been the most effective tool to appraise inflation because the CMB 
is the earliest electromagnetic "snapshot" of the universe, and as such, probes the 
universe in a particularly pristine state - before gravitational and electromagnetic 
processing. Since gravity is the weakest of the four fundamental forces, gravitational 
radiation (the inflationary GWB) probes much farther back - to '--^ 10~^^ sec after 
the Big Bang. 

CMB polarization is generated by both scalar and tensor perturbations. The 
two are related in all models of inflation since both are generated by fluctuations of 
the same quantum field. Inflation may have occurred at energies too low to det ect 
(or not at all) but the scalar-tensor relation provides a powerful consistency checlP^ 
- insurance against a false-positive claim. 

All inflationary observables are determined by the inflaton potential. The pro- 
cess is invertible, affording the opportunity to reconstruct, or significantly constrain, 
the inflaton potential from CMB measurements of the GWB. Assuming the infla- 
tionary paradigm is correct, this provides a window into physics at energy scales 
below 10^^ GeV even if future CMB observations produce a null result. Without 
observational co nstraint s, the number of potential inflationary models has swelled 
to daunting levelJ^^I^. Upcoming experiments should be able to rule out many of 
these models, even with a null result as we discuss in section [TUl 
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CMB temperature anisotropy measurements can determine the amplitude of the 
primordial power spectrum, As, and the scalar power-spectral index, n^, far better 
than CMB polarimeters but only polarization sensitive experiments can measure 
the key-prediction of inflation - the tensor-to-scalar ratio, r, if it is below r ~ 0.3, 
which is near the current 2a upper limits as discussed in ^ A joint detection of Ag, 
Us, and r will serve to reconstruct the inflaton potential (as well as completely rule 
out the ekpyrotic and variable-speed-of- light models). This observation requires an 
understanding of the effects of gravitational waves on the CMB, to which we now 
turn. 



4. Primordial Gravitational Wave CMB Perturbations 

Gravitational waves perturb the metric tensor describing the geometry of the early 
universe. The general, perturbed, metric of a flat FLRW universe can be written as 



ds^ = (77) [-di]'^ + {Sij + h^j)dx^dx^] . (1) 

The components htj (77, x) represent the gravitational waves and can be ex- 
panded into spatial Fourier harmonics e^*'' '', 

h., iv, x) = / E (1^) [^^ ^'"'^ +hl (77) e-'^ - C] , (2) 

where k is a time-independent wave vector, and k — {Sijk^yy^'^ . The wave 
number k defines the wavelength measured in units of laboratory standards by 

s 

A = 2iTa/k. The polarization tensors Pij (k) have different forms depending on 
whether the y re present gravitational waves, rotational perturbations, or density 
perturbation^^. In the case of gravitational waves the polarization tensors can be 
expressed in terms of two mutually orthogonal unit-vectors (1, m) lying in the plane 
of the wave front (i.e. perpendicular to k), 

1 2 

Pij= hlj — 'mirrij, Pij= kvij + Ijirii, (3) 
and obey the conditions 

kj 5'^ = 0, k, = 0, 
P,, (k) p' (k) = 26s' s- (4) 

For a classical gravitational wave field, the quantities and in Eq. |21 are 
arbitrary complex (conjugate) numbers. 

The Fourier expansion allows us to reduce the problem of evolution of the 
perturbed gravitational field to the evolution of mode functions hk{r]) for each 
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individual mode k. For a single Fourier component the amplitude of the primordial 
gravitational wave obeys the following wave equatiori'*^: 



drf a drj drj ' 

This equation ignores the damping effects of anisotropic stress provided by 
cosmological neutrinos^. For fcry <C 1 (wavelength larger than the cosmological 
horizon) the last term in the above equation can be dropped out, and the amplitude 
of the gravitational wave is "frozen" [hk ~ const). For krj ^ I (wavelength smaller 
than the horizon) the solutions for hk are damped plane waves. Since larger k 
perturbations enter the horizon earlier, their contribution to the tensor spectrum is 
more heavily damped. 

Equation [S] allows us to study the gravitational wave perturbations in two 
cosmological epochs: 1) in the matter dominated epoch, governed by matter with 
effective equation of state p = 0, where the scale factor behaves as a{r]) cx r]"^ and 
2)in the radiation dominated epoch, governed by the effective equation of state 
p = e/3, with a{r]) cx rj. 

In the general case with cosmic scale factor a(r/), Eq.O might not allow for an 
analytical solutioi i. F or a mixed matter-radiation Universe, the scale factor is given 
by the expression 1^ 



a{ri) = 



41 



H 



1 + + 



?7(?7 + 2?7*), 



(6) 



where rj^ = (\/2 — l)'7egi a-nd ijeq is the time of matter-radiation equality. The 
participating constants have been chosen such that the value of the scale factor 
at the present epoch rjo is a(?7o) = '21 h, which leads to the value of the time at 



present epoch 7/0 



/2. With such a convention, the wave 



(l-2ry,) + v/l + 4r;2 
whose wavelength, A, today is equal to present Hubble radius, carries the con- 
stant wavenumber kn = 47r. The scale factor has two asymptotic regions: rj <^ rj^f 
corresponding to the radiation dominated epoch with a{r]) oc 77, and rj ^ corre- 
sponding to the matter dominated epoch with a(r/) cx 77^. 

An elegant method to analyze the evolution of gravitational waves in a such 
a universe is to smoothly approximate the above scale factor with power-law scale 
factors, for which we can analytically solve the equation (jS)). The scale factor (Eq. 
is well approximated by 



a{ri) = 21h ■ {rj ~ rj^f , V** < V- 



(7) 



Imposing continuity of 0(77) and a'(77) at matter-radiation equality, 77**, fully deter- 
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mines all participating constants. 

_ 1 



The general solution to Eq. |S1 for the scale factor (Eq. [71) is given by 

hk{^) = <^ (9) 

The constants Cfc and are determined by the evolution of the gravitational 
waves before radiation domination. Since the radiation dominated era is believed 
to have been preceded by inflation the values of the coefficients Ck and Dk, which 
give the spectral characteristics of the gravitational wave fie ld, are determined by 
the physics of inflation and initial conditions (see for example^ I I ). It follows that 
due to the initial stage of rapid expansion Ck « —Dk- 

The fact that Ck ~ —Dk shows that the gravitational wave modes, k, are 
(almost) standing waves at the radiation epoch. To find the coeSicients Ak and 
Bk in the solution hkiv) and dhk{ri) / drj must be joined continuously at the 
transition point rj = rj^^ Is^I 

The above analytical approximation @ sufflciently well approximates the so- 
lution to the equation (0) in the case of the scale factor given by ®. Numeri- 
cal calculations show that the above analytical approximations work very well for 
wavenumbers k satisfying kArj^q < 1, where Arjeq is the characteri stic time scale of 
change from the radiation-dominated era to matter-dominated er 

5. CMB Polarization Observables 

Having considered the underlying cosmological behavior of tensor perturbations we 
now turn to the imprint of tensor gravitational waves on CMB polarization using the 
equation of radiative transport. To begin we consider a polarized electromagnetic 
wave with angular frequency, co: 

E = Eyo sin{ujt — Sy)y + E^o sin(wi — 5x)y^- 

The polarization state of is characterized by the Stokes parameters: /, Q, U, and V: 

I = Iy+I^, 

with ly = (Eyf^) and Ix — {EIq). I is the total intensity of the radiation, and is 
always positive. The other Stokes parameters are defined as 

Q = ly — Ix and 

U = 2EyoExo cos{6y - Sx) 

V = 2Ey(,Exo sin{Sy - bx) 
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where Q and U quantify the hnear polarization of the wave, and V quantifies the 

degree of circular polarization. The polarization fraction is 11 = , and 

the polarized intensity is /poi = 11 x /. The Stokes parameters comprise a symbolic 
vector i introduced by Chandrasekhail^ and related to the Stokes parameters in 
the following way: 

ly 
U 

\vj 

Thomson scattering only produces linear CMB polarization, implying V — 0, 
so we will only consider the symbolic 3-vector: I = 

\u. 

Polarized radiation in the presence of cosmological metric perturbations is rep- 
resented as state vector describing the occupation numbers of polarized radiation^ 2 











no + n^Sn ~ uq 




;■) 


+ 5h 











(10) 



The Boltzmann equation of radiative transfer written in terms of n(r/, x" , 6, ip) 



is: 



dn 
dr] 



dn 



dn dv 

dv drj 



g(n-J) 



and 



47r./o JO 



27r _ 



P{e, V, d', ^p')n{ri, , ^, 61', V') sin O'dO'diP', 



(11) 



(12) 



where q = axNed, a is the cosmological scale factor, J is the "scattering" or "coUi- 
sional" term which is a function of the angular variables only, with primed angular 
variables corresponding to the photon direction before scattering and unprimcd 
angular variable s co rresponding to the photon direction after scattering. P is the 
scattering matrijJ^, ax is the Thomson cross section, and A^e is the comoving num- 
ber density of free electrons. The coupling of the gravitational waves to the radiation 
is manifested in the first term on the right side of Eq. ^2 

V dr] 2 drj 
('\ 

Let us introduce fio = no 1 corresponding to unpolarized isotropic ther- 

Vo/ 

mal radiation (zero-th order approximation), with riQ = [exphv/kbT — 1]~^, which 
depends only on the photon frequency v and corresponds to the Planck spectrum. 
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The angular symmetry of P requires: 

J(no) = — / / P(6',V,6'',V')nosin6l'de''dV' = no. (13) 

It should be mentioned that the form of the scattering integral in H13|l assumes that 
the chosen reference frame is comoving with the scattering electrons. In the case of 
scalar perturbations an additional Doppler term mavarise due to the movement of 
electrons with respect to the chosen reference frame^. In the case of gravitational 
waves this Doppler term does not arise, since it is always possible to choose our 
synchronous reference frame to be comoving with the scattering electrons. 

In Eg. 1101 (5n is the first order correction to the uniform, isotropic, and unpolar- 
ized radiation described by ng. This perturbation is comprised of an (unpolarized) 
term due to angular anisotropy of the photon distribution, iiA, and a polarized 
term, fin, so (5n = ha + fin- 

The anisotropic and polarized components are functions of conformal time, ry, 
co-moving spatial coordinates, a;", photon frequency, v, and photon propagation 
direction specified by the unit vector e°'{6,(f)) with polar angle, 9, and azimuthal 
angle, (p. All of the polarization-specific cosmological phenomena discussed in this 
review can be described by their effect on 6n. These effects, if measurable, can be 
used to evaluate and refine the standard cosmological model. Since CMB polar- 
ization depends on anisotropy, it probes all the same underlying physics and, in 
addition, several cosmological effects are only observable via CMB polarization and 
not CMB anisotropy. For this reason, CMB polarization is a valuable cosmological 
tool. The penalty we pay, however, is that the Thomson scattering which produces 
CMB polarization has a fairly inefficient coupling to cosmological perturbations. We 
will evaluate the cosmological phenomena to which CMB polarization is sensitive 
in subsequent sections. 

We retain only the zero and first order perturbation terms in ha/s- Since ^ 
is of the first order, we can replace ^ by in Eq. ^2 (i/g is the unperturbed 
frequency). This implies that the frequency dependence of both the polarization 
and anisotropy is given by the same factor!^ 

^ vq dnp 
no dvo 

We now concentrate on the first-order terms. In the following we will identify 
J(n) = Ji((5fi). An arbitrary gravitational wave can be considered as a linear su- 
perposition of plane gravitational waves. Due to the linear nature of the problem 
the anisotropy and polarization generated by an arbitrary gravitational wave is the 
linear superposition of anisotropy and polarization generated by plane gravitational 
waves. After spatial Fourier transformation, the first-order Boltzmann equation be- 
comes: 



— ikfii,Sn{T], k) = - 2^~a^ 



g(r;)[M(7?,k)-J(,7,k)]. (14) 
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where, fXk = j^" , k — |k| and 0k is the azhiiuthal angle of e" in the plane per- 
pendicular to the vector k. In the case of a single gravitational wave we can choose 
our spherical coordinate system in such a way that cos 9 — fij^ and (j) — 4>k (we shall 
omit the index k in /i and (j) when such an omission does not lead to confusion). 



6. CMB polarization due to gravitational waves. An analytical 
approach. 

We now focus on the solutions to Eq. E| and with our attention restricted to 
primordial gravitational waves (tensor perturbations). The process of generation of 
polarization occurs after radiation-matter equality. In this case the source in Eq. 1141 
has the following form: 



^T^e-e'^ = ^7(1 " m') cos 2cj^S{v, k), (15) 

where S{ri, k) = '^'^^j^^^ is the gravitational wave source term. 

For a plane gravitational wave perturbation with wavevector k, the synrbolic 
vector (5n(77,k), describing anisotropy and polarization, can be presented aS^ 



Y{ri,fi,k){l - ^2) 



cos 20 + P{ri, ^, k) 



(1 + ^2) cos 20 
-(1 ^2) cos 20 
4/i sin 20 



(16) 



For krj 1, a and (3 do not depend on fj.^. Equation 1161 will allow for recon- 
struction of the polarization and temperature power spectra, which can be derived 
from maps of the CMB's polarization. 

Substituting H15|l and (|16|l into the Boltzmann equation, H14|l . we obtain the 
following system of coupled integro-differential equations for a and (3 



Kv. M, k) + {q- ikfi)(3{rj, fi, k) = —q{vi)I{r], k), 



(17) 



iiv, k) + {q- ikfi)^{r], fi, k) = S{'q, k), 
where ^(rj, /i, k) = a{r], /i, k) + l3{rj, /i, k) and 



/(77, k) 



(1 + /i'')'/3(r?,/i', A:) - -(1 - /i'2)2^(ry,M', fc) 



(18) 



(19) 



Here 



q{r]) ^ axNeaXeiri) = 



ar^BPcXeir]) 



where ax is the Thomson optical depth, f^s is the baryon fraction, pc is the critical 
density, nip is the proton mass, a is the cosmological scale factor, and Xe{ri) is the 
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ion ization fraction. For Xe(ri) during decoupling we use the Peebles fitting function 
El. 

The formal solution of Eq. El is 

/3(r?, k) = e"('')+*'='^'' P fc)e^*^-^"da;, (20) 

^0 

where T(77)is optical depth: 

r{v) - r qiv')dv'- (21) 

and 



V 



<^{Ti,k) = ^I{Ti,k)q{in)e'-^^\ (22) 
16 

The function ^r],k{v) depends primarily on the epoch and duration of decou- 
pling, and is not sensitive to the details of T{rij^. 
Taking into account that from Eq. E| 

^{r,, k) = e"('')+''=^'' /" 5(x, A:)e-"(")-''=''"dx, (23) 



(24) 



Jo 

we obtain for $ the following integral equation: 



3 r 

$(77, k) = $o(?7, k) + —q{r]) / $(a;, k)K+{'q - x, k)da 



where 



and 



Mv,k)^~^q{il) I dxS{x,k)e-^^-^K^{rj-x,k), (25) 



K±{x,k) = J ^dfi{l±^ 



fi ) cos kfix. (26) 

We remind the reader that our goal is to obtain the present value of the symbolic 
vector 6n, from Eq. El To do so requires that $(77, /c) be found for every k. The 
vector Sh is found by by using equations (jJOJ, 123|l and substituting into Eq. I|16|l . 
setting rj = TjQ and T(r/o) = 0. Once the symbolic vector <5n is found, we proceed 
to determine the multipole expansion of the GWB-induced polarization, which will 
eventually be used to construct statistical estimators of observable quantities such 
as the temperature and polarization power spectra. 

7. Multipole expansion of Anisotropy and Polarization due to 
gravitational waves 

This section is largely based on Seljak & Zaldarriag The components of the 
symbolic vector n(/i, 0) are related to the fundamental polarization tensor Patip-, 4') 
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in the following way 



1 

2h^ 



U 



U /- 



ni 



■a:b 



(28) 
(29) 
(30) 



(27) 

From the polarization tensor Pab{fJ-,(t>) we can form three independent scalar 
fields corresponding to anisotropy and E, B modes of polarization in the following 
way : 

E (/^, 0) = [gab (m, 0) n (m, 0) - 2Pa6 (M, 

i? iii, <t>) = e"^, • [gab (m, 0) • n (M, - 2Pa6 (/i, (t>)t"' : 

where repeated indices imply summation, and 

/(l-M^)-i 

gab ifJ-, 0) = 

V (i-M^), 

is the 2-metric on a unit sphere in coordinates {fi,(j)), and covariant derivatives (;) 
are effected with respect to this metric. 

In terms of our variables a (/i, k) and /3 (/i, k) (see the above definitions of 
T, and B (|28I30^ give 

T (/i, fc) = 7 (l - /i^) a (^, fc) • cos 2(j), 



£;(M,0,fc) = -7(l-M') 

B (/i>, fc) = 7 (1 - m') 



(3 (^, fc) cos2(/), 
/3 {fi, k) sin 20, 



t— 8 — 

The multipole expansion coefficients for anisotropy, and E, B modes of polar- 
ization, in terms of the introduced scalars T,E and i?, are defined as 



Ytm (m, 0) • T (a*, (j}, k) 




(Yl,^{^i,4,) ) ■E{fi,<i>,k) 



(31) 
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where Y/^m (/i, 4>) are the ordinary spherical harmonics. From H31() it foUows that the 
multipole coefficients vanish for m ^ ±2. Taking into account the form expressions 
for a and (3 from H20|l and (j2Sl, the above expressions can be rewritten as: 







'([+2)! m 

(1-2)1 C2 
1(1-1) 



^2 J jiiO - -^ji-iiC) 



aF,mik) = -{-i)'l{52,ra - 5-2,m) V^(2Z + 1) j dj] Hfc (ry) 



where 



(l-l) 



C = k(vo - v), 



and 



Tfc(77) = 5fc(r/)e-^('')-$fe (77), 
and ji(C) are the spherical bessel functions. 

8. Anisotropy and polarization generated by a random field of 
gravitational waves 

The GWB can be considered as an isotropic random superposition of plane waves 
with the following correlation relationship: 



(CfcCk') = (ckCk') ^6ss'6(k-k), (ckCj^, 



= 



(32) 



where (...) means averaging over realizations. Below we assume power law spectrum 
for Cosmological Gravitational Wave Background, i.e. 



4TrP\hk(v)\' 



i,nT — 1 



with Tlx — 1 corresponding to flat (scale-invariant) Zel'dovich-Harrison spectrum. 

Let us consider two scalar fields X(jjL,(j)) and X(iJ,,(j)) (where X and X is any 
pair of the scalars T(ii,(f)), E(iJL,(f) and B(ii,(j))). Then their cross correlation is 
defined as 

T^^(Ho) = ^ / dfidn'd(l)d(l)' X 



x5 (^0 - A^m' - V(1-m')(1-M'2)cos (0 - 0')) {X(^i, 0)A(Ai', ^')). 

(33) 
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Following Fourier decomposition of cf)) 



s=l,2 



34) 



(27r)3/2 

and the correlation relations H32|) , the correlation function is presentable in the form 
r-^"^(^o) = 47r / k^dk r^^{no,k), (35) 



where T-^-^ {iiq, k) is the correlation function for a single wave 



1 



r^^(/io,fc) d^i V d0 del)' 



6 (a^o - M/^' - V(1-m')(1-A*'')cos (0 - 0')) (m', 

(36) 

The correlation function being a function of a single angle 6*0 = cos~^ /zq, can 
be expanded over Legendre polynomials to give the following correlation function 
and power spectrum 



r^*(Mo,fc)-EQ^^(fc)Pi(/io) 



/=0 

1 



21 + 1 



(37) 



(38) 



Using the addition theorem for Legendre polynomials 



Pi 



(1 - m2)(1 - A*'2)cos(0-0') = Pz(m)P,(a.')+E 



m=0 



(/-m)! 
{l + m)\ 



Pr(M)^r(M)cosm(0-0), 



we get 



1 1 

2/ + 14^ 



d/io Pli^J'o)d^ d^' dcj) dip' 



/ d^i dii' dcp dcp' Pi ( ^Ifl' + yir^TI^Kw^^cos (</) - 0')) 

X(pi, 0, k)X*{n', (/)', k), 



21 



1 

TjTY ^ ^l,rn{k) ■ Xi'„^{k), 



m— — l 



where 



Xl^rn{k) ^ f dfl (Yi^yn{^l,(f)))* ■ X(pL,cf),k) , 



(39) 



(40) 
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(41) 



Thus we get for the power spectrum from a single gravitational wave 




1 



Ex X 



,x* 

•Lm 



(42) 



2/ + 1 



The angular power spectrum from a superposition of gravitational waves is 
calculated from the above expression by integrating over all the wave numbers k 



The power spectra is most conveniently determined as coefficients of a Legendre 
polynomial expansion with index I rather than an expansion in Fourier modes with 
wavenumber k. In a flat universe a perturbation of comoving wavelength k~^ at the 
comoving distance of the last scattering surface (LSS) subtends an angle 9 ~ 1/k. 
On the other hand £ ^ 1/9 ior small 6. This simple consideration shows that the 
main contribution to Cf^ with a given £ comes from k ^ £. 

Figure El shows Cf^ calculated using the CAMB^"^- code, which is based on 
CMBFASlI^, compared to our analytical predictions. The peak multipole of the 
spectrum £peak — 90 is robust to changes in the inflationary dynamics. The am- 
plitude of the B-mode power spectrum, is determined by r alone, and its spatial 
structure is determined solely by the age of the universe at last scattering. This 
single-parameter dependence makes the B-mode polarization the most robust probe 
of inflation. While scalar, or mass-energy, perturbations are amplified by gravity, the 
tensor-GWB is not. Direct detection today (redshift z = 0) by, e.g., LIGO/LISA is 
essentially impossible since, like the CMB, the energy-density of gravitational waves 
dilutes (redshifts) for waves inside the horizon at, or before, last-scattering as the 
universe expands. However the GWB imprints curl- mode polarization on the CMB 
at the surface of last-scattering (z ~ 1100). Therefore, the GWB energy-density 
from sub-horizon scale waves at last-scattering was at least one trillion times larger 
than it is now, which motivates the use of the last-scattering surface as the perhaps 
the best "detector" of the primordial, inflationary-generated, GWB 

Our analytic calculations based on the solution of the Boltzmann equation 
provide clear insight into the underlying physics of polarization and are in good 
agreement with the resuhs of CMBFAST and CAMB. 

9. Current CMB Polarization Results 

Figures Q and |21 show current detections of the E-mode (grad-mode) polar iza- 
tion and the polarization-temperature cross correlation, (TE). Currently, DASl^, 
BOOMERANd^, CbJ^, and CAPMAI^Sni have detected E-mode polarization, 
however there are no current detections of the E-mode polarization for £ < 100 




(43) 
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0.0 



00 



Fig. 1 . Measurements of the gradient or E-mode polarization power spectrum Cf ■ The sohd line 



is the polarization power spectrum for the WMAP best fit cosmological model, with r = 0.17 



EH 



where the signature of gravitational waves will be manifest. Only DASI and 
BOOMERANG have detected the grad-mode polarization at 100 < £ < 400. 

Both reionization and gravitational waves imprint the polarization of the CMB 
at large scales. The primary effects of reionization are encoded in the grad-mode and 
temperature polarization cross-correlation at £ < 50. In the absence of detections, 
the most stringent constraint in this range of multipoles is currently Cf < 8 ij.K 
at 95% confidence reported by POLABf^for 2 < < 20 (assuming no B-modes). 
WMAP^ reports a large number of highly significant detections, especially at low- 
i due to WMAP's ability to map the full sky. CBI, DASI, and BOOMERANG 
have also detected the cross-correlation spectrum, mainly at smaller angular scales 
than WMAP. A complete description of reionization will require detections of CMB 
E-mode polarization for £ < 50. An ancillary benefit of reionization is that it 
boosts the primary curl-mode power spectrum significantly near £ — 10. Due to 
reionization, a more stringent limit on the tensor-to-scalar ratio, r, in the presence 
of lensing can be obtained than that calculated in Ref. I61L 1621 
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Fig. 2. Measurements of the temperature-polarization cross-correlation power spectrum ■ 
Th e solid line is the power spectrum for the WMAP best fit cosmological model, with t = 0.17 

EH 



The tensor to scalar ratio, r, is defined here as 

_ Aliko) 



A|(fco) 



(44) 



where A|j(A;o) and A^(A;o) are the amplitudes of the scalar and primordial power 
spectra, evaluated at some pivot wavenumber feg . There are other ways of defining 
the tensor to scalar ratio that appear in the literature. A second definition, T/S, 
is defined as the ratio between the tensor and scalar contributions to the temper- 
ature anisotropy's quadrupole component T/S — C2_tensl^2 seal- relationship 
between r and T / S depends in a complicated way on the cosmology. For the pa- 
rameter set used in this paper, T/S « 0.5r where the pivot wavenumber used is 
0.05 Mpc-i. 

CMB temperature anisotropy a^one can only detect the tensor to scalar ratio 
if r > 0.3 due to cosmic variance^. The limit from WMAP data alone is roughly 
four-times larger than this and approximately two-times larger with the inclusion of 



February 5, 2008 12:11 WSPC/INSTRUCTION 

keating'polnarev3march2006'FINAL 



FILE 



18 Brian G. Keating, Alexander G. Polnarev, Nathan Miller, and Deepak Baskaran 

external datcP^. Thus we are close to the ultimate constraints achievable by CMB 
temperature anisotropy observations alone. 

Further progress requires detection of the B-mode signal. Currently, there are no 
detections of the B-mode polarization as shown in figureOl POLAR, CBI, DASI, and 
BOOMERANG all provide upper limits to the B-mode polarization. The current 
generation of CMB polarimeters should be able to provide detect or provide much 
better upper limits to the curl-mode polarization within the next several years. 




10 100 1000 

/ 

Fig. 3. Measurements of the curl-mode polarization power spectrum Cf^ . Th e solid line is the 
power spectrum for the WMAP best fit cosmological model with t = O.ITI^^ and r = 0.1. The 
dotted line are the analytical results of this paper, and the solid curve which peaks at £ ~ 1000, 
is the B-mode spectrum produced by large scale structure lensing the primary CMB E-mode 
polarization. 

The gravitational lensing of the E-mode polarization into B-mode polarization 
provides a source of contamination to the primordial B-mode signal. For a noise free, 
full sky experiment in which the lensing is treated as noise, R ef. ICTl found that this 
contamination sets a detectability limit of rum > 10-4|611d21 j g jf f]^^ energy scale 
of inflation is larger than 3 x lO^^GeV^. Measurements of 21-cm radiation may be able 
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to provide a way to de-lens the curl-mode measurements, providing an detectable 
inflationary energy scale limit lower than one using CMB-only measurements. This 
energy scale may be as low as 3 x 10^^ GeV^. Work by Seljak & Hirata (2004) 
1^ indicate that the primary lensing signal can be removed to a level that makes 
r = 10"^ detectable, allowing detection of inflationary GWB at energy scales < 10^^ 
GeV using the CMB only. 

10. Constraining Inflation with Upcoming CMB Observations 

To illustrate the power of upcoming observations to constrain inflationary parame- 
ters we consider a hypothetical polarimeter with 1° resolution and a system sensitiv- 
ity {NETsys) = 70/iKsi/2 using the COSMOMC^^: software package. The system 
sensitivity for polarization is \/2 times highe, NEQ = \/2NET. These detector re- 
quirements are well within the reach of the current generation of CMB polarimeters. 
By observing 2.4% of the sky, this experiment could feasibly detect the curl mode 
signal created by the gravitational wave background (at t < 100) for r as low as 
r = 0.12 at 95% confidence with no priors. This represents an order of magnitude 
improvement over the WMAP only results of r = 1.28 at 95% confidence, with 
no priors. We plot the results of our simulations along with predictions for several 
classes of inflationary models in figure 01 The predictions of the various models of 
inflation were generated using the method described below. 

Detecting the tensor-to-scalar ratio, r, will allow the possibility of distinguishing 
between several classes of inflation, such as negative curvature models, small positive 
curvature models, etc. As an example we consider the toy experiment described 
earlier. As mentioned, it can detect r with 95% at 0.12. The ris—r contours compared 
to the different classes of inflation models are plotted in Figure 01 assuming the toy 
experiment detects r — 0.12 and Ug = 0.98. The contours for the toy experiment 
were calculated using COSMOMC. 

The predictions for the different classes of inflation were generated following 
the model in Ref. l67[l^ and the different inflation models were separated according 
to Ref. The observables r and Ug were evaluated at some specific e-folding, N, 
of inflation, not a speciflc wavenumber, k. The relation between N and k requires 
a detailed model of reheating, which has some uncertainty. This uncertainty is 
marginalized over by calculating the observables at an e-fold randomly drawn from 
40 to 70. The inflationary flow equations were truncated at sixth order and the 
observables were calculated to second order in slow roll. Each class of model has 
a unique color in Fig. 0] as described in the caption, and each realization of that 
model is indicated by a point in the r — Ug plane. 

11. Conclusions 

We have developed an analytic method to generate the predictions of the imprint 
of gravitational waves on the CMB. Using phenomenological models of inflation we 
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Fig. 4. Shown are ID marginalized probabilities and 2D joint likelihood contours for Us and r 
from the toy experiment described in the text. The purple dots are negative curvature models. The 
red dots are small positive curvature models. The green dots are intermediate positive curvature 
models. The black dots are large positive curvature models. WMAP has set an upper limit on r of 
r < 1.28 with no priors. The toy experiment is capable of detecting r = 0.12 with 95% confidence 
with no priors— more than an order of magnitude below the WMAP limits. 

predict both the CMB polarization spectra and the derived inflationary parame- 
ters: the tcnsor-to-scalar ratio and spectral index of the scalar perturbations. The 
combination of the later two observables allows for reconstruction of the dynamics 
of inflation. With these predictions in hand we have shown that upcoming CMB 
polarization observations will be able to detect or constrain the cosmological GWB 
and hence, inflation itself. These new technological advances now position obser- 
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vational cosmology at the threshold of an exhilarating era - one in which CMB 
polarization data will winnow down inflation's vast model-space and test models of 
the early universe at energy scales approaching the GUT-scale; nearly one trillion 
times higher energy than accessible from particle accelerators. 
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